Directed, Elliptic and Triangular Flows in Asymmetric Heavy Ion 

Collisions 



M. Bleicher 1 ' 3 , K. A. Bugaev 2 , P. Rau 13 , A. S. Sorin 4 , J. Steinheimer 1 ' 3 , and 

H. Stocker 1 - 5 



1 Frankfurt Institute for Advanced Studies (FIAS), 60438 Frankfurt, Germany 
2 Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 03680 

Kiev, Ukraine 

,_i 3 Institut fiir Theoretische Physik, Goethe-Universitat, 60438 Frankfurt, Germany 

4 Joint Institute for Nuclear Research (JINR), 141980 Dubna, Moscow Region, Russia 
5 GSI Helmholtzzentrum fur Schwerionenforschung, 64291 Darmstadt, Germany 

I— s 

emails: bleicher@th.physik.uni-frankfurt.de bugaev@th.physik.uni-frankfurt.de 
rau@th.physik.uni-frankfurt.de sorin@theor.jinr.ru steinheimer@th.physik.uni-frankfurt.de 

h.stoecker@gsi.de 

-t— > 
I 

Ch Abstract 

In this paper we propose to thoroughly investigate asymmetric nuclear collisions both in 
J> the fixed target mode at the laboratory energy below 5 GeV per nucleon and in the collider 

mode with a center of mass energy below 11 GeV per nucleon. Using the UrQMD transport 

^ model, we demonstrate a strong enhancement of directed and elliptic flow coefficients for the 

midcentral asymmetric nuclear collisions compared to symmetric collisions. We argue that 

\q such an enhancement is due to the disappearance of the nuclear shadowing effect on the side 

of the smaller projectile nucleus. An analysis of the energy and centrality dependencies of 
the directed, elliptic and triangular flows at midrapidity shows us their sensitivity to the 
details of the employed model of hadronic interaction. In general, the flow patters found 

• — for asymmetric nuclear collisions have a very rich and complicated structure of energy and 

centrality dependencies compared to the flows found for symmetric collisions and are worth 
to be investigated experimentally. The directed, elliptic, and triangular flow coefficients are 
computed for target nuclei containing high density fluctuations and thoroughly compared 
with that ones obtained in the absence of such fluctuations. 
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1. Introduction. The experimental study of strongly interacting matter has reached a 
decisive moment: it is hoped that the low energy heavy ion collisions programs performed 
at the CERN SPS and the BNL RHIC and two new programs that are planned to begin in 
a few years at NICA (JINR, Dubna) and FAIR (GSI, Darmstadt) will allow the heavy ion 
community to locate the mixed phase of the deconfinement phase transition and to discover a 
possible (tri)critical endpoint. The major experimental information [1, 2, 3, 4, 5] is provided 



by the measurements of particle yields, one particle momentum spectra, two particle correla- 
tions, and the Fourier components of the collective hadronic flow [6] known as ui-coefncient 
(directed flow), ^-coefficient (elliptic flow), and t>3-coefficient (triangular flow). Although 
the great success of experiments at the BNL RHIC [1, 2] and at the CERN LHC [3, 4, 5] 
proved the high efficiency of modern experimental methods, their results also clearly demon- 
strated that the heavy ion collisions programs at RHIC, SPS, NICA and FAIR energy range 
are not simpler and they require further development of both new experimental approaches 
and far more sophisticated theoretical models in order to reach their goals. Therefore, in 
view of new opportunities opening with the Nuclotron program Baryonic Matter@Nuclotron 
(BMN) [7, 8] which will start at JINR (Dubna) in 2014, we would like to address here some 
new physical issues that can be studied at laboratory energies of 2 — 5 AGeV for a wide range 
of colliding nuclei in the fixed target mode. In future, they can also be further investigated 
at the accelerators of new generation like NICA (JINR, Dubna) and FAIR (GSI, Darmstadt) 
at a center of mass energy up to 11 AGeV. 

This range of energies was thoroughly investigated in the past at the GSI SIS and the 
BNL AGS experiments, but only for symmetric nuclear collisions A+A. In this paper we 
demonstrate by using the Ultrarelativistic Quantum Molecular Dynamics transport model 
UrQMD [9, 10] that one of the major tasks of the low energy programs at JINR and GSI could 
be a systematic study of directed, elliptic, and triangular flows for non-central asymmetric 
nuclear collisions (ANC), i.e. for non-central A + B reactions with l<B<C7l (see Fig. 1). 
This aspect of heavy ion physics was not yet systematically explored. Until today there 
were only a few works reported for ANC [11, 12, 13]. Moreover, these works were completed 
before a systematic investigation of the Fourier components [6] of hadronic flow was proposed. 
Therefore, here we argue that at lab. energies of about 2—10 AGeV a systematic study of 
vi, V2 and vs Fourier-coefficients of the azimuthal particle distributions measured in non- 
central ANC with a special choice of impact parameter may help to essentially improve our 
understanding of the hadronic matter equation of state at high densities. In addition, under 
certain conditions, highest baryonic charge densities can be reached in these experiments [14], 
what could help clarifying the question, whether the mixed quark-hadron phase [7, 8, 15, 16] 
or the predicted chiral quarkyonic phase [17, 18, 19, 16] can be formed in this energy range. 

First theoretical predictions in this energy region were made very recently [14] and they 
indicate a very rich picture of the physical phenomena in ANC to be investigated, e.g. the 
formation of high density Mach shock waves and their impact on the spatial distribution of 
the produced particles and heavier clusters. The main purpose of this work is to formulate 
the most promising physical issues for the ANC of the BMN program at JINR and to work 
out the strategy of their experimental exploration. 

This work is organized as follows. In the next section we discuss the difference of the 
directed, elliptic, and triangular flows obtained in ANC and in symmetric nuclear collisions. 
The third section is devoted to the analysis of the directed, elliptic and triangular flows that 
may develop, if high density fluctuations occur inside the target nucleus. The last section 
contains our conclusions. 

2. The unusual properties of ANC. ANC have some history since they were sug- 
gested long ago (see, for instance, [11, 20, 21, 22, 23]), but in those days the analysis of 
Fourier coefficients of the azimuthal particle distributions was not even suggested. Since 
the experiments on collisions of identical nuclei is simpler, ANC were forgotten for awhile. 
After the first work [6] on the analysis of the Fourier coefficients v\ , V2 , • • • was published 
this method has become a powerful tool for experimental studies of the evolution process of 
symmetric heavy ion collisions at high energies [24, 25, 26, 27, 28], i.e. of two identical nuclei 
(A + A). Thus, the quark scaling of the i>2 dependence on the transverse momentum pt [28] 
clearly demonstrated the partonic source of elliptic flow of hadrons at RHIC energies, while 
the triangular flow is reflecting the correlations that appear at the early stage of collisions 
[29]. 
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Fig. 1. Schematic picture of an asymmetric nuclear collision of nuclei A (large circle) with B (small circle) 
shown in the transverse plane. The common area of two nuclei is shown for very peripheral collisions 
(panel a)), for semi-peripheral collisions (panels b) and c)) and for the most central collision (panel d)). 
Clearly, for semi-peripheral collisions the effect of the shadowing of particle motion to the right hand side 
of the nucleus B is very weak. 

A principally new element of the ANC compared to symmetric nuclear collisions is the 
generation of strong and asymmetric gradients of energy density and baryonic density at 
the initial stage of the collision process along with the stronger flow from the target in the 
direction of the projectile nucleus for a specific choice of the impact parameter values. In 
contrast to symmetric collisions, in ANC the reflection symmetry between the left nucleus A 
and the right nucleus B (see panels a)-c) in Fig. 1) is broken from the very beginning. This 
leads to an entirely different shape and location of the overlap region between the colliding 
nuclei which, in its turn, results in different flow patterns compared to symmetric A + A 
collisions. Indeed, if the impact parameter value is close to bANC ~ 1.1 fm • (Az — Bs ) ± 1 fm 
(see the panels b) and c) of Fig. 1) and if the size of the target nucleus is chosen close to -Bs « 
\Ai , then there is sufficient room to vary the impact parameter in the experiment and to 
select values which are close to bANC using an event-by-event analysis. In this case ANC allow 
us to scan the interior of the target nucleus, as well as to study in detail the variety of surface 
phenomena such as the surface formation of light nuclear fragments like deuterons, tritons, 
and heliums nuclei, the emission of high pj- pions [30], the isospin dependence of hadron 
surface emission, conical emission due to Mach shock wave formation in the target [14], and 
so on. 

Consider the change in the directed flow first. Although the approximate dependence 
of the so-called flow parameter F [31] on Az + E>3, i.e. on the collision time, was found 
in the symmetric nuclear collision experiments [32], however, a similar dependence of the 
directed flow in the ANC is much less certain [26] and has to be studied. In addition, in 
symmetric nuclear collisions with a lab. energy between 1 and 10 AGeV the effect of nuclear 
shadowing plays an important role and causes the negative value of pionic x>\ [25] . Since both 
the peripheral and the most central ANC (see panels a) and d) of Fig. 1) are very similar to 
the corresponding symmetric nuclear collisions due to a similar geometry, one can expect a 
similar behavior of their v\ coefficients, i.e. v^ s ~ vf , while for the semi-peripheral ANC (see 
panels b) and c) in Fig. 1) one expects a different situation. This expectation is supported 
by the experimental analysis of [12] of pionic flow directed to the in-plane OX-axis and the 
one directed oppositely. Therefore, the semi-peripheral ANC would allow one to essentially 
reduce the effect of shadowing on one side of the formed fireball what could help to finally 
clarify the source of the strong antiflow of pions [30, 33] with p? < 500 MeV observed at 
SIS energies [34] and predicted earlier [35, 36]. Moreover, one can hope to study a rich 
structure of shock waves in hadronic matter (from conical emission to viscosity) occurring 
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Fig. 2. Energy and centrality dependence of the vf s coefficient of charged particles from Ne+Au 
collisions as calculated with the UrQMD model with potentials (left panel) and without them (right 
panel). For more details see the text. 
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Fig. 3. Energy and centrality dependence of the vf s coefficient found by the UrQMD model with 
potentials for all charged particles (left panel) and for positive pions (right panel) for the ANC Ne+Au. 



during semi-peripheral ANC [14]. 

In order to demonstrate the new possibilities of the ANC we performed the analysis of 
Vi S using the UrQMD model. This model is one of the most successful transport models 
and was supplemented by in- medium potentials [37, 38, 39] in order to reproduce the full 
complexity of the flow patterns in the low energy range. For this work, we use the UrQMD 
model developed in [9, 10, 37, 38, 39] which was thoroughly tested on the available data in a 
wide range of collision energies. Fig. 2 shows the energy and centrality dependence of the vf s 
coefficient for all charged particles from the ANC Ne+Au at longitudinal midrapidity in the 
equal velocity frame of the colliding nuclei. Such a system is convenient to compare the results 
with both the data obtained in symmetric nuclear collisions and with UrQMD simulations 
reported earlier [39]. As one can see from Fig. 2 the vf s coefficient of charged particles is 
essentially non-zero at \y\ < 0.1 for ANC, whereas it is zero for y = in symmetric nuclear 
collision since it is an odd function of the center of mass rapidity. Also Fig. 2 demonstrates 
that that inclusion of the potentials in the UrQMD model [37, 38, 39] is important since their 
absence leads to roughly double decrease of the directed flow coefficient. A comparison of the 
charged particles directed flow with that one of positive pions shows a complete similarity 
in the energy and centrality dependence (see Fig. 3). The behavior of negative pions is 
also similar. Figs. 2 and 3 demonstrate a gradual decrease of the directed flow (left-right 
asymmetry) with the increase of the collision energy. 

Let us now consider the qualitative change of the elliptic flow (in-plane) coefficient be- 
havior in ANC. The experimental data show that at lab. energies of about 3.3 AGeV the 




elliptic flow in symmetric collisions is zero, i.e. i>2 = 0) an d its value slowly grows with the 
colliding energy [27]. The reason for those small values of the elliptic flow around this en- 
ergy is that the slow moving remnants of the target (projectile) nucleus A (B) shadow the 
motion of particles from the collision zone directed from nucleus A to nucleus B (see Fig. 1). 
The change of the sign at 3.3 AGeV signifies the transition from the out-of-plane to the 
in-plane elliptic flow which was predicted long ago [40]. 

Due to a similar geometry both the peripheral and the most central ANC are alike the 
corresponding symmetric nuclear collisions, and therefore, in this case one has to expect 
almost identical behavior of their V2 coefficients, i.e. v^ s ~ • On the other hand, the semi- 
peripheral collisions (panels b) and c) in Fig. 1) produce an entirely different situation: in 
this case, if the smaller nucleus B hits the larger nucleus A at its boundary, not too far above 
or below its surface, i.e for impact parameter values close to bANC, then the whole situation 
concerning the shadowing is changed. At these impact parameter values there is almost no 
shadowing on one side and, hence, one can expect that in those ANC the v£ S coefficient 
can essentially be enhanced at lab. energies between 2 and 5 AGeV. In addition, one can 
expect that some number of slow moving particles with an initial momentum directed from 
nucleus B to nucleus A is reflected backwards and creates an additional flow in the direction 
of the no-shadowing side of the collision region leading to positive values of ■ These 
new features of the ANC elliptic flow can be seen in Fig. 4. For Eigb = 2 AGeV both the 
charged particles and pions have positive v^ s coefficient for semi-peripheral collisions, while 
for central and peripheral collisions the pionic v^ s is negative. As the collision energy is 
increased, the maximum of the f^ 5 coefficient gradually moves to larger impact parameter 
values and gets less pronounced. 

In this work we report the first simulations of the triangular flow coefficient for the 
low collision energy range. Recent analyses [29, 1, 2, 3, 4, 5, 41] clearly demonstrated an 
importance of the t>3 coefficient for an elucidation of the experimental information and here 
we confirm this fact. Although the v^ s signal found for E\. A ^ = 2 — 10 AGeV is much weaker 
than that one determined at RHIC energies, the structure of triangular flow at low collision 
energies is definitely much richer as one can see from Fig. 5. Indeed, for E\ a b = 2 — 3 AGeV 
the charge particles (negative pions) predict that the v$ s coefficient has a maximum in semi- 
peripheral collisions whereas v$ s is negative at impact parameter values of b ~ 1 fm (b < 1 
fm) and b > 6 fm. For E^ = 4 — 5 AGeV the maximum of the v% s coefficient is still clearly 
seen, although it gets wider, whereas for larger collision energies this maximum is not that 
clearly visible especially for pions (see the right panel in Fig. 5). 

3. The dense spot effect on flow patterns. The existence of high density fluctuations 
in ordinary nuclei is debated for a long time [42, 43] along with their possible applications in 
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Fig. 5. Energy and centrality dependence of the v% s coefficient found by the UrQMD model with 
potentials for charged particles (left panel) and for negative pions (right panel) in the ANC Ne+Au. 
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Fig. 6. Same as in Fig. 3, but for the dense spot configuration (see text for details). 



the context of the heavy ion reactions [44, 45]. Here we would like to recapitulate this issue 
and to study the influence of the high density fluctuations on the vf s , v^ S and v% s coefficients 
for ANC. To estimate this effect, we randomly put 20 nucleons from the Au-nucleus in a 
narrow Gaussian distribution inside of the ordinary Monte Carlo sampled nucleus. This 
gives us a fluctuating dense spot initial configuration for the target nucleus. Such a mild 
assumption is far from the extremely high densities discussed with respect to flucton [42, 45]. 
The dense spot is fixed in the center of the Au target with an offset in x-direction of half 
of the impact parameter between the two colliding nuclei. The condition that there is an 
overlap between the dense spot of the radius 1.1 B 3 fm and the projectile nucleus of the same 
radius is 6 — 1.1 1?3 fm < | + 1.1 Bz fm. That means, for impact parameters of the Ne+Au 

reaction of b < AABz fm ~ 12 fm such an overlap always exists. 

The results for the v^ s , v^ S and v^ s coefficients are presented in Figs. 6-8. Comparing 
Figs. 3 and 6, one sees that the directed flow for collisions with the target exhibiting the 
dense spot and for collisions without it have almost an identical centrality dependence for 
the energies £"i a b = 2 — 4 AGeV. However, for > 5 AGeV the v ± s coefficient with the 
dense spot is just 50 % of the vf s coefficient without dense spot. Thus, there is a non- 
monotonic behavior of the directed flow coefficient with respect to the collision energy. From 
the analysis of energy density we conclude that it indicates the change of regimes from the 
dominance of the target break up process at E^b = 2 — 4 AGeV to a strong compression 
and more intense thermalization of the reaction zone for £"i a b > 5 AGeV. Qualitatively 
a similar picture is valid for the elliptic flow coefficient of the dense spot configuration: 
for Si a b = 2 — 4 AGeV the coefficients with the dense spot are similar to that ones 
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Fig. 7. Same as in Fig. 4, but for the dense spot configuration 
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without it, whereas for > 5 AGeV the elliptic flow coefficient with the dense spot has 
different centrality dependence (compare Figs. 4 and 7). The latter means that, in contrast 
to the case of the dense spot absence, for > 5 AGeV the elliptic flow coefficient of the 
charged particles is close to zero, i.e. ~ 0, for b < 3 fm while its peak is now located 
at b = 6 fm. On the other hand, the coefficient of positive pions exhibits a W-shape 
for E\ab = 2 AGeV, whereas for E^ > 4 AGeV its centrality dependence is similar to the 
elliptic flow coefficient of protons and charged particles. Also, one has to admit that for 
Eiab > 4 AGeV the maximal amplitudes of the elliptic flow of both the charged particles, 
protons and pions for the dense spot configurations (Fig. 7) are not reduced in strength 
compared to the ANC without dense spot (Fig. 4). Thus, the 50 % reduction of the directed 
flow for the dense spot configurations at E^ > 4 does not lead to a dramatic change of the 
corresponding elliptic flow at these energies. 

Again the essential change in the centrality dependence of the vf s coefficient for the 
dense spot configuration is seen at E^ < 4 AGeV (see Fig. 8): the clear maximum located 
at about b = 3 — 4 fm for the no dense spot case (see Fig. 5) disappears in the dense spot 
configurations. Only for energies E\gb > 8 AGeV one can see some increase of the maximal 
value of v% s coefficient compared to the no dense spot case. 

4. Conclusions. Using the UrQMD model including potentials [37, 38, 39] here we 
perform a comprehensive analysis of the directed, elliptic and triangular flow coefficients for 
asymmetric nuclear collisions (ANC) for the first time and suggest to explore this promising 
aspect of heavy ion physics in the energy range of the BMN program at JINR and the FAIR 
heavy ion program. We argue that a specific choice of the impact parameter should lead to 
the disappearance of the nuclear shadowing effect on the side of smaller colliding nucleus and 
show that such an effect leads to an essential enhancement of the directed and elliptic flows 



of hadrons in ANC compared to the symmetric nuclear collisions. 

Our analysis shows that the centrality and energy dependencies of the v ± s , v^ s and v% s 
coefficients found for ANC are richer and more complicated compared to symmetric nuclear 
collisions. Furthermore, we find that these flow patterns are very sensitive to the details of 
the employed interaction which can be used both for fine tuning of transport codes and for 
elucidation of the essential features of the hadron interaction in the medium. 

Also, we demonstrate that the possible existence of high density fluctuations in the large 
target nucleus may be verified using ANC, since they allow one to scan the target interior 
by a smaller projectile and to draw physical conclusions from the centrality and the energy 
dependencies of directed, elliptic and triangular flows on an event-by-event basis. The found 
non-monotonic energy dependence of the directed flow indicates a higher degree of compres- 
sion reached in ANC with the dense spot configuration at lab. energies above 5 AGeV. 

All these new results allow us to hope that the low energy ANC will soon become a power- 
ful tool of theoretical and experimental studies of the strongly interacting matter properties. 
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